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C/5 Abstract 

The suppression and collective flow of open and hidden heavy flavor in ultra-relativistic heavy ion collisions is studied with the 
partonic transport model Boltzmann Approach to MultiParton Scatterings (BAMPS). Charm and bottom quarks interact with the 
rest of the medium via elastic scatterings with an improved Debye screening and a running coupling. Heavy flavor electron data 
from RHIC can only be described if the cross section is multiplied with K = 3.5 to mimic the contribution from radiative processes. 
With the same value for K we compare to elliptic flow data of D mesons at the LHC and find a good agreement. However, the 
nuclear modification factor of D mesons, non-prompt J/if/, muons, and electrons at the LHC is underestimated in our calculations. 
In addition, we present first BAMPS results on the suppression and elliptic flow of J /if/ mesons, which are treated in the same 
framework. 
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1. Introduction 

Experimental measurements at the Relativistic Heavy Ion Collider (RHIC) Q] ,2] and Large Hadron Collider 
(LHC) |3 1 reveal that in ultra-relativistic heavy ion collision a hot and dense medium is created which exhibits exciting 
properties like collective flow and jet quenching. This medium consists of deconfined quarks and gluons and, hence, 
is called quark gluon plasma (QGP). Heavy quarks in particular are an interesting probe to study the properties of 
^vq the medium since they can only be produced in initial hard parton scatterings or in the early phase of the medium 

t— I evolution due to their large mass [4]. Therefore, they traverse the medium right from the beginning for a rather long 

time, interact with other medium particles, participate in the flow and lose energy. 
• <-h Due to flavor conservation the number of heavy quarks stays constant during hadronization. Open heavy flavor 

mesons can then be measured (directly or indirectly), which reveals important information about the heavy quark 
distribution and hence the properties of the QGR At the LHC, ALICE has reconstructed D mesons directly and 
measured their nuclear modification factor Raa [5 1 and elliptic flow V2 @. In addition, also the B meson Raa has been 
measured via non-prompt J /if/ by the CMS collaboration [7]. More indirect are the measurements of heavy flavor 
muons [8] and electrons [9, 1 ] which have contributions from both D and B mesons. All measurements show a large 
energy loss of heavy quarks in the QGP and a strong participation in the collective flow. 

Complimentary to open heavy flavor mesons with only one heavy quark are hidden heavy flavor mesons like J '/if/ 
which consist of a charm and an anti-charm quark. Lattice calculations indicate that J /if/ states can survive in the 
QGP to some extent [11]. The most interesting observable here is the nuclear modification factor of J /if/. However, in 
contrast to the open heavy flavor Raa it is not dominated by energy loss, but by the total number of J /if/ in the medium. 
Nevertheless, this observable is also strongly influenced by cold nuclear matter effects like shadowing which makes it 
a challenging task for theoretical explanations. 
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Figure 1 . Elliptic flow vi (left) and nuclear modification factor Raa (right) of D mesons, non-prompt J I if/, muons, and electrons as a function of 
transverse momentum for Pb+Pb collisions at LHC together with data 1 5 6 10 24] . As a note, the centrality classes of the Raa data are not always 
the same. The impact parameter of the BAMPS calculations is matched to the corresponding centrality class to be able to compare to the data. 



2. The partonic transport model BAMPS 

Boltzmann Approach to MultiParton Scatterings (BAMPS) ifTZl [T3l is a partonic transport model which solves 
the Boltzmann equation for on-shell particles and pQCD interactions. For light partons all 2 — > 2 as well as 2 <-> 3 
processes are implemented. On the heavy flavor sector so far only binary scatterings of heavy quarks with light 
partons are active, namely, gQ — > gQ and qQ — > qQ. The divergent t channel of those processes is regularized with a 
screening mass p which is determined by matching energy loss calculations with the Born cross section to results from 
hard thermal loop calculations lfT4lfl6l . The comparison of both results shows that the screening mass p is of the order 
of the Debye mass m D , p 2 = Kin 2 D with k = l/(2e) * 0.2. Furthermore, the running of the coupling is explicitly taken 
into account. More details on both can be found in Ref. [ 16 1. We note that currently for light partons neither a running 
coupling nor an improved Debye screening is employed, which we plan, however, to do in the future. Inelastic heavy 
flavor processes ifTTll are currently being studied in BAMPS and results will be presented soon. 

In addition, J/tfr mesons have been added to BAMPS which can be dissociated by a gluon, J/iff + g — > c + c [ 18 1, 
and regenerated via the back reaction of this process. Furthermore, cold nuclear matter effects, that is in BAMPS 
shadowing, nuclear absorption, and the Cronin effect, are introduced, which is outlined in detail in Ref. |[T9ll . 



3. Results on open heavy flavor 

The running coupling and the improved screening procedure, which reproduces the energy loss from hard thermal 
loop calculations, effectively enhances the heavy quark cross section with the medium. Quantitative comparisons 1 16, 
I20H231 show that elastic processes are responsible for a large fraction of the energy loss of heavy quarks. However, 
they are not able to reproduce the data of the nuclear modification factor or elliptic flow of any heavy flavor particle 
species. This is not too surprising since we expect that radiative 2 — > 3 processes also play an important role and 
that both processes together should account for the measured suppression and flow. The quantitative contribution of 
radiative processes will be studied in an forthcoming investigation. In this paper we mimic the radiative influence 
by effectively increasing the elastic cross section by a factor K — 3.5 which fits the heavy flavor elliptic flow data 
from RHIC 11231 . Simultaneously, the Raa data at RHIC can also be described with the same parameter. Whether 
this effective description is valid and radiative contributions really boil down to simply multiplying the elastic cross 
section by a constant factor will be studied with BAMPS in the future. 

Figure [T] gives an overview of the BAMPS results on the elliptic flow and nuclear modification factor of all heavy 
flavor particles which can be measured at the LHC. The experimental data is also plotted wherever it is available. 
The BAMPS predictions for the elliptic flow of D mesons are in good agreement with the data points over the whole 
transverse momentum range. This hints that the effective binary cross section with K - 3.5, which fits the RHIC 
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Figure 2. Left panel: Raa of J /if/ at mid-rapidity |v| < 0.35 and forward rapidity 1.2 < y < 2.2 for Au+Au collisions at RHIC as a function of the 
number of participants, together with experimental data 1 25 1. Right panel: Elliptic flow of J/if/ with data 1 26 ]. For comparison the charm quark vj 
is also shown. 



data, can also describe the elliptic flow at LHC, although the experimental errors are still too large to draw any 
definite conclusion. The v-i of D mesons is considerably larger than that of non-prompt J /if/ due to the mass difference 
of charm and bottom quarks. The elliptic flow of electrons and muons stemming from both D and B mesons is 
somewhere in between. 

The nuclear modification factor of all heavy flavor species is at the lower edge of the error bars or even below it. 
Therefore, the suppression in BAMPS with K — 3.5 is too strong, although it agrees with the RHIC Raa of heavy 
flavor electrons. However, at RHIC the error bars especially at high transverse momentum are rather large. Therefore 
it could be possible that also at RHIC scaled elastic collisions alone cannot describe both vi and Raa simultaneously 
and radiative collisions must be taken into account. 



4. Suppression and elliptic flow of J /iff 

First BAMPS results on the suppression of J/ift mesons at RHIC were presented in Ref. |[T9l . In Fig.[2]the Raa 
of J/if/ at RHIC is depicted as a function of the number of participants. For this calculation a formation time of to = 
0.6 fm is employed for J /if/ mesons to prevent early melting when the temperature cannot be properly defined. The 
experimental data of the Raa at mid-rapidity is well described. However, our results at forward rapidity underestimate 
the suppression for central and semi -central events. 

In this section we want to focus primarily on the elliptic flow of J/if/. Experimental measurements at RHIC showed 
that the V2 of J /if/ is very small |26 1. This is in contradiction with the regeneration picture where the flow of the charm 
quarks should be transferred to the J /iff. BAMPS is an ideal framework to study this in more detail since it reproduces 
the D meson flow (cf. previous section) and also allows recombination of charm quarks to J /if/. The right panel of 
Fig.|2]shows that in BAMPS even with regeneration the elliptic flow of all J /if/ is compatible with the data. 

To estimate the effect of regeneration on flow in the following qualitative study, we consider only the elliptic flow 
of secondary J/if/, which are created in the medium from the recombination of two charm quarks. Initial J /if/ has 
(nearly) vanishing V2 and must be taken into account if one wants to compare to data. In Fig. [3] we show the number of 
secondary J fi/i as a function of time and their elliptic flow as a function of transverse momentum for two scenarios. In 
one scenario the J /if/ melting temperature is large. Hence, secondary J /if/ can be formed already in the early phase of 
the medium evolution where the temperature is still large. In the other scenario we choose a small melting temperature 
which only allows secondary J /if/ production at late times. 

J / if/ from the first scenario only have a small V2 since they are produced early and their constituent charm quarks 
could not build up a sizeable flow. The second scenario corresponds more to the coalescence picture. The J/if/ V2 
is shifted to twice as large values at twice as large pj compared to the charm flow, v^ipj) m 2v c 2 (pj/2). Here the 
fraction of initial J /if/ at large pj must be considerably larger than secondary J /if/ to be compatible with the data. 
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Figure 3. Time evolution of the number of secondary J /if/ (normalized) and their elliptic flow as a function of transverse momentum at the end of 
the QGP phase at RHIC for small and large J/ifi melting temperature Tj. For comparison, the charm vi and experimental data on J /if/ [26] are also 
shown. As a note, the contribution of initial J /if/ with vanishing \>2 is not depicted. 



5. Conclusions 

We presented recent BAMPS results on the elliptic flow and nuclear modification factor of several heavy flavor 
particles, namely D mesons, non-prompt J /if/, muons, and electrons. Using only elastic collisions of heavy quarks 
with medium particles with an improved Debye screening inspired by hard thermal loop calculations and a running 
coupling the RHIC heavy flavor electron data can only be reproduced with scaling the binary cross section with 
K = 3.5. With the same parameter we find a good agreement with the D meson elliptic flow data at LHC. However, 
the Raa of D mesons, non-prompt J /if/, and muons are underestimated. Furthermore, we showed some preliminary 
results on the suppression and elliptic flow of J /i/r with a focus on flow of secondary J /if/ that stem from charm quarks 
which develop flow within the same framework. 
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